A functional gene module identification algorithm in gene expression data based on genetic algorithm and gene ontology.

Since genes do not function individually, the gene module is considered an important tool for interpreting gene expression profiles. In order to consider both functional similarity and expression similarity in module identification, GMIGAGO, a functional Gene Module Identification algorithm based on Genetic Algorithm and Gene Ontology, was proposed in this work. GMIGAGO is an overlapping gene module identification algorithm, which mainly includes two stages: In the first stage (initial identification of gene modules), Improved Partitioning Around Medoids Based on Genetic Algorithm (PAM-GA) is used for the initial clustering on gene expression profiling, and traditional gene co-expression modules can be obtained. Only similarity of expression levels is considered at this stage. In the second stage (optimization of functional similarity within gene modules), Genetic Algorithm for Functional Similarity Optimization (FSO-GA) is used to optimize gene modules based on gene ontology, and functional similarity within gene modules can be improved. Without loss of generality, we compared GMIGAGO with state-of-the-art gene module identification methods on six gene expression datasets, and GMIGAGO identified the gene modules with the highest functional similarity (much higher than state-of-the-art algorithms). GMIGAGO was applied in BRCA, THCA, HNSC, COVID-19, Stem, and Radiation datasets, and it identified some interesting modules which performed important biological functions. The hub genes in these modules could be used as potential targets for diseases or radiation protection. In summary, GMIGAGO has excellent performance in mining molecular mechanisms, and it can also identify potential biomarkers for individual precision therapy.

MVIP: multi-omics portal of viral infection.

Virus infections are huge threats to living organisms and cause many diseases, such as COVID-19 caused by SARS-CoV-2, which has led to millions of deaths. To develop effective strategies to control viral infection, we need to understand its molecular events in host cells. Virus related functional genomic datasets are growing rapidly, however, an integrative platform for systematically investigating host responses to viruses is missing. Here, we developed a user-friendly multi-omics portal of viral infection named as MVIP (https://mvip.whu.edu.cn/). We manually collected available high-throughput sequencing data under viral infection, and unified their detailed metadata including virus, host species, infection time, assay, and target, etc. We processed multi-layered omics data of more than 4900 viral infected samples from 77 viruses and 33 host species with standard pipelines, including RNA-seq, ChIP-seq, and CLIP-seq, etc. In addition, we integrated these genome-wide signals into customized genome browsers, and developed multiple dynamic charts to exhibit the information, such as time-course dynamic and differential gene expression profiles, alternative splicing changes and enriched GO/KEGG terms. Furthermore, we implemented several tools for efficiently mining the virus-host interactions by virus, host and genes. MVIP would help users to retrieve large-scale functional information and promote the understanding of virus-host interactions.

Analysis of bronchoalveolar lavage fluid metatranscriptomes among patients with COVID-19 disease.

To better understand the potential relationship between COVID-19 disease and hologenome microbial community dynamics and functional profiles, we conducted a multivariate taxonomic and functional microbiome comparison of publicly available human bronchoalveolar lavage fluid (BALF) metatranscriptome samples amongst COVID-19 (n = 32), community acquired pneumonia (CAP) (n = 25), and uninfected samples (n = 29). We then performed a stratified analysis based on mortality amongst the COVID-19 cohort with known outcomes of deceased (n = 10) versus survived (n = 15). Our overarching hypothesis was that there are detectable and functionally significant relationships between BALF microbial metatranscriptomes and the severity of COVID-19 disease onset and progression. We observed 34 functionally discriminant gene ontology (GO) terms in COVID-19 disease compared to the CAP and uninfected cohorts, and 21 GO terms functionally discriminant to COVID-19 mortality (q < 0.05). GO terms enriched in the COVID-19 disease cohort included hydrolase activity, and significant GO terms under the parental terms of biological regulation, viral process, and interspecies interaction between organisms. Notable GO terms associated with COVID-19 mortality included nucleobase-containing compound biosynthetic process, organonitrogen compound catabolic process, pyrimidine-containing compound biosynthetic process, and DNA recombination, RNA binding, magnesium and zinc ion binding, oxidoreductase activity, and endopeptidase activity. A Dirichlet multinomial mixtures clustering analysis resulted in a best model fit using three distinct clusters that were significantly associated with COVID-19 disease and mortality. We additionally observed discriminant taxonomic differences associated with COVID-19 disease and mortality in the genus Sphingomonas, belonging to the Sphingomonadacae family, Variovorax, belonging to the Comamonadaceae family, and in the class Bacteroidia, belonging to the order Bacteroidales. To our knowledge, this is the first study to evaluate significant differences in taxonomic and functional signatures between BALF metatranscriptomes from COVID-19, CAP, and uninfected cohorts, as well as associating these taxa and microbial gene functions with COVID-19 mortality. Collectively, while this data does not speak to causality nor directionality of the association, it does demonstrate a significant relationship between the human microbiome and COVID-19. The results from this study have rendered testable hypotheses that warrant further investigation to better understand the causality and directionality of host-microbiome-pathogen interactions.

Weighted Gene Co-Expression Network Analysis to Identify Potential Biological Processes and Key Genes in COVID-19-Related Stroke.

The purpose of this research was to explore the underlying biological processes causing coronavirus disease 2019- (COVID-19-) related stroke. The Gene Expression Omnibus (GEO) database was utilized to obtain four COVID-19 datasets and two stroke datasets. Thereafter, we identified key modules via weighted gene co-expression network analysis, following which COVID-19- and stroke-related crucial modules were crossed to identify the common genes of COVID-19-related stroke. The common genes were intersected with the stroke-related hub genes screened via Cytoscape software to discover the critical genes associated with COVID-19-related stroke. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis for common genes associated with COVID-19-related stroke, and the Reactome database was used to annotate and visualize the pathways involved in the key genes. Two COVID-19-related crucial modules and one stroke-related crucial module were identified. Subsequently, the top five genes were screened as hub genes after visualizing the genes of stroke-related critical module using Cytoscape. By intersecting the COVID-19- and stroke-related crucial modules, 28 common genes for COVID-19-related stroke were identified. ITGA2B and ITGB3 have been further identified as crucial genes of COVID-19-related stroke. Functional enrichment analysis indicated that both ITGA2B and ITGB3 were involved in integrin signaling and the response to elevated platelet cytosolic Ca2+, thus regulating platelet activation, extracellular matrix- (ECM-) receptor interaction, the PI3K-Akt signaling pathway, and hematopoietic cell lineage. Therefore, platelet activation, ECM-receptor interaction, PI3K-Akt signaling pathway, and hematopoietic cell lineage may represent the potential biological processes associated with COVID-19-related stroke, and ITGA2B and ITGB3 may be potential intervention targets for COVID-19-related stroke.

Screening of Gene Expression Markers for Corona Virus Disease 2019 Through Boruta_MCFS Feature Selection.

Since the first report of SARS-CoV-2 virus in Wuhan, China in December 2019, a global outbreak of Corona Virus Disease 2019 (COVID-19) pandemic has been aroused. In the prevention of this disease, accurate diagnosis of COVID-19 is the center of the problem. However, due to the limitation of detection technology, the test results are impossible to be totally free from pseudo-positive or -negative. Improving the precision of the test results asks for the identification of more biomarkers for COVID-19. On the basis of the expression data of COVID-19 positive and negative samples, we first screened the feature genes through ReliefF, minimal-redundancy-maximum-relevancy, and Boruta_MCFS methods. Thereafter, 36 optimal feature genes were selected through incremental feature selection method based on the random forest classifier, and the enriched biological functions and signaling pathways were revealed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Also, protein-protein interaction network analysis was performed on these feature genes, and the enriched biological functions and signaling pathways of main submodules were analyzed. In addition, whether these 36 feature genes could effectively distinguish positive samples from the negative ones was verified by dimensionality reduction analysis. According to the results, we inferred that the 36 feature genes selected via Boruta_MCFS could be deemed as biomarkers in COVID-19.

Mechanism of COVID-19-Related Proteins in Spinal Tuberculosis: Immune Dysregulation.

Purpose: The purpose of this article was to investigate the mechanism of immune dysregulation of COVID-19-related proteins in spinal tuberculosis (STB). Methods: Clinical data were collected to construct a nomogram model. C-index, calibration curve, ROC curve, and DCA curve were used to assess the predictive ability and accuracy of the model. Additionally, 10 intervertebral disc samples were collected for protein identification. Bioinformatics was used to analyze differentially expressed proteins (DEPs), including immune cells analysis, Gene Ontology (GO) and KEGG pathway enrichment analysis, and protein-protein interaction networks (PPI). Results: The nomogram predicted risk of STB ranging from 0.01 to 0.994. The C-index and AUC in the training set were 0.872 and 0.862, respectively. The results in the external validation set were consistent with the training set. Immune cells scores indicated that B cells naive in STB tissues were significantly lower than non-TB spinal tissues. Hub proteins were calculated by Degree, Closeness, and MCC methods. The main KEGG pathway included Coronavirus disease-COVID-19. There were 9 key proteins in the intersection of COVID-19-related proteins and hub proteins. There was a negative correlation between B cells naive and RPL19. COVID-19-related proteins were associated with immune genes. Conclusion: Lymphocytes were predictive factors for the diagnosis of STB. Immune cells showed low expression in STB. Nine COVID-19-related proteins were involved in STB mechanisms. These nine key proteins may suppress the immune mechanism of STB by regulating the expression of immune genes.

Discovering Common Pathophysiological Processes between COVID-19 and Cystic Fibrosis by Differential Gene Expression Pattern Analysis.

Coronaviruses are a family of viruses that infect mammals and birds. Coronaviruses cause infections of the respiratory system in humans, which can be minor or fatal. A comparative transcriptomic analysis has been performed to establish essential profiles of the gene expression of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) linked to cystic fibrosis (CF). Transcriptomic studies have been carried out in relation to SARS-CoV-2 since a number of people have been diagnosed with CF. The recognition of differentially expressed genes demonstrated 8 concordant genes shared between the SARS-CoV-2 and CF. Extensive gene ontology analysis and the discovery of pathway enrichment demonstrated SARS-CoV-2 response to CF. The gene ontological terms and pathway enrichment mechanisms derived from this research may affect the production of successful drugs, especially for the people with the following disorder. Identification of TF-miRNA association network reveals the interconnection between TF genes and miRNAs, which may be effective to reveal the other influenced disease that occurs for SARS-CoV-2 to CF. The enrichment of pathways reveals SARS-CoV-2-associated CF mostly engaged with the type of innate immune system, Toll-like receptor signaling pathway, pantothenate and CoA biosynthesis, allograft rejection, graft-versus-host disease, intestinal immune network for IgA production, mineral absorption, autoimmune thyroid disease, legionellosis, viral myocarditis, inflammatory bowel disease (IBD), etc. The drug compound identification demonstrates that the drug targets of IMIQUIMOD and raloxifene are the most significant with the significant hub DEGs.

A Comparative Transcriptomic Analysis Reveals That HSP90AB1 Is Involved in the Immune and Inflammatory Responses to Porcine Deltacoronavirus Infection.

PDCoV is an emerging enteropathogenic coronavirus that mainly causes acute diarrhea in piglets, seriously affecting pig breeding industries worldwide. To date, the molecular mechanisms of PDCoV-induced immune and inflammatory responses or host responses in LLC-PK cells in vitro are not well understood. HSP90 plays important roles in various viral infections. In this study, HSP90AB1 knockout cells (HSP90AB1KO) were constructed and a comparative transcriptomic analysis between PDCoV-infected HSP90AB1WT and HSP90AB1KO cells was conducted using RNA sequencing to explore the effect of HSP90AB1 on PDCoV infection. A total of 1295 and 3746 differentially expressed genes (DEGs) were identified in PDCoV-infected HSP90AB1WT and HSP90AB1KO cells, respectively. Moreover, most of the significantly enriched pathways were related to immune and inflammatory response-associated pathways upon PDCoV infection. The DEGs enriched in NF-κB pathways were specifically detected in HSP90AB1WT cells, and NF-κB inhibitors JSH-23, SC75741 and QNZ treatment reduced PDCoV infection. Further research revealed most cytokines associated with immune and inflammatory responses were upregulated during PDCoV infection. Knockout of HSP90AB1 altered the upregulated levels of some cytokines. Taken together, our findings provide new insights into the host response to PDCoV infection from the transcriptome perspective, which will contribute to illustrating the molecular basis of the interaction between PDCoV and HSP90AB1.

Augmentation of Transcriptomic Data for Improved Classification of Patients with Respiratory Diseases of Viral Origin.

To better understand the molecular basis of respiratory diseases of viral origin, high-throughput gene-expression data are frequently taken by means of DNA microarray or RNA-seq technology. Such data can also be useful to classify infected individuals by molecular signatures in the form of machine-learning models with genes as predictor variables. Early diagnosis of patients by molecular signatures could also contribute to better treatments. An approach that has rarely been considered for machine-learning models in the context of transcriptomics is data augmentation. For other data types it has been shown that augmentation can improve classification accuracy and prevent overfitting. Here, we compare three strategies for data augmentation of DNA microarray and RNA-seq data from two selected studies on respiratory diseases of viral origin. The first study involves samples of patients with either viral or bacterial origin of the respiratory disease, the second study involves patients with either SARS-CoV-2 or another respiratory virus as disease origin. Specifically, we reanalyze these public datasets to study whether patient classification by transcriptomic signatures can be improved when adding artificial data for training of the machine-learning models. Our comparison reveals that augmentation of transcriptomic data can improve the classification accuracy and that fewer genes are necessary as explanatory variables in the final models. We also report genes from our signatures that overlap with signatures presented in the original publications of our example data. Due to strict selection criteria, the molecular role of these genes in the context of respiratory infectious diseases is underlined.

Bioinformatics/network topology analysis of acupuncture in the treatment of COVID-19: response to methodological issues.

Coronavirus disease 2019 (COVID-19) is a global pandemic and there is an urgent need to discover the therapy for COVID-19. In our original article, we first obtained the target proteins of acupuncture and related target genes of COVID-19 by searching English and Chinese databases, then Gene Ontology biological processes and enrichment analysis were performed on the overlapping targets of acupuncture with COVID-19. Moreover, the compound-target and compound-disease-target network was constructed. This is an innovative attempt to predict the potential benefits of acupuncture treatment for COVID-19. In this letter, we answered reader Zheng's comments.

New Approach for Targeted Treatment of Mild COVID-19 by Honeysuckle through Network Pharmacology Analysis.

OBJECTIVE: To investigate the potential pharmacological value of extracts from honeysuckle on patients with mild coronavirus disease 2019 (COVID-19) infection. METHODS: The active components and targets of honeysuckle were screened by Traditional Chinese Medicine Database and Analysis Platform (TCMSP). SwissADME and pkCSM databases predict pharmacokinetics of ingredients. The Gene Expression Omnibus (GEO) database collected transcriptome data for mild COVID-19. Data quality control, differentially expressed gene (DEG) identification, enrichment analysis, and correlation analysis were implemented by R toolkit. CIBERSORT evaluated the infiltration of 22 immune cells. RESULTS: The seven active ingredients of honeysuckle had good oral absorption and medicinal properties. Both the active ingredient targets of honeysuckle and differentially expressed genes of mild COVID-19 were significantly enriched in immune signaling pathways. There were five overlapping immunosignature genes, among which RELA and MAP3K7 expressions were statistically significant (P < 0.05). Finally, immune cell infiltration and correlation analysis showed that RELA, MAP3K7, and natural killer (NK) cell are with highly positive correlation and highly negatively correlated with hematopoietic stem cells. CONCLUSION: Our analysis suggested that honeysuckle extract had a safe and effective protective effect against mild COVID-19 by regulating a complex molecular network. The main mechanism was related to the proportion of infiltration between NK cells and hematopoietic stem cells.

Integrative multiomics and in silico analysis revealed the role of ARHGEF1 and its screened antagonist in mild and severe COVID-19 patients.

COVID-19 is a sneaking deadly disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The rapid increase in the number of infected patients worldwide enhances the exigency for medicines. However, precise therapeutic drugs are not available for COVID-19; thus, exhaustive research is critically required to unscramble the pathogenic tools and probable therapeutic targets for the development of effective therapy. This study utilizes a chemogenomics strategy, including computational tools for the identification of viral-associated differentially expressed genes (DEGs), and molecular docking of potential chemical compounds available in antiviral, anticancer, and natural product-based libraries against these DEGs. We scrutinized the messenger RNA expression profile of SARS-CoV-2 patients, publicly available on the National Center for Biotechnology Information-Gene Expression Omnibus database, stratified them into different groups based on the severity of infection, superseded by identification of overlapping mild and severe infectious (MSI)-DEGs. The profoundly expressed MSI-DEGs were then subjected to trait-linked weighted co-expression network construction and hub module detection. The hub module MSI-DEGs were then exposed to enrichment (gene ontology + pathway) and protein-protein interaction network analyses where Rho guanine nucleotide exchange factor 1 (ARHGEF1) gene conjectured in all groups and could be a probable target of therapy. Finally, we used the molecular docking and molecular dynamics method to identify inherent hits against the ARHGEF1 gene from antiviral, anticancer, and natural product-based libraries. Although the study has an identified significant association of the ARHGEF1 gene in COVID19; and probable compounds targeting it, using in silico methods, these targets need to be validated by both in vitro and in vivo methods to effectively determine their therapeutic efficacy against the devastating virus.

Identification of the Key Immune-Related Genes in Chronic Obstructive Pulmonary Disease Based on Immune Infiltration Analysis.

Purpose: Chronic obstructive pulmonary disease (COPD) is a major cause of death and morbidity worldwide. A better understanding of new biomarkers for COPD patients and their complex mechanisms in the progression of COPD are needed. Methods: An algorithm was conducted to reveal the proportions of 22 subsets of immune cells in COPD samples. Differentially expressed immune-related genes (DE-IRGs) were obtained based on the differentially expressed genes (DEGs) of the GSE57148 dataset, and 1509 immune-related genes (IRGs) were downloaded from the ImmPort database. Functional enrichment analyses of DE-IRGs were conducted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses and Ingenuity Pathway Analysis (IPA). We defined the DE-IRGs that had correlations with immune cells as hub genes. The potential interactions among the hub genes were explored by a protein-protein interaction (PPI) network. Results: The CIBERSORT results showed that lung tissue of COPD patients contained a greater number of resting NK cells, activated dendritic cells, and neutrophils than normal samples. However, the fractions of follicular helper T cells and resting dendritic cells were relatively lower. Thirty-eight DE-IRGs were obtained for further analysis. Functional enrichment analysis revealed that these DE-IRGs were significantly enriched in several immune-related biological processes and pathways. Notably, we also observed that DE-IRGs were associated with the coronavirus disease COVID-19 in the progression of COPD. After correlation analysis, six DE-IRGs associated with immune cells were considered hub genes, including AHNAK, SLIT2 TNFRRSF10C, CXCR1, CXCR2, and FCGR3B. Conclusion: In the present study, we investigated immune-related genes as novel diagnostic biomarkers and explored the potential mechanism for COPD based on CIBERSORT analysis, providing a new understanding for COPD treatment.

Differential Co-Expression Network Analysis Reveals Key Hub-High Traffic Genes as Potential Therapeutic Targets for COVID-19 Pandemic.

Background: The recent emergence of COVID-19, rapid worldwide spread, and incomplete knowledge of molecular mechanisms underlying SARS-CoV-2 infection have limited development of therapeutic strategies. Our objective was to systematically investigate molecular regulatory mechanisms of COVID-19, using a combination of high throughput RNA-sequencing-based transcriptomics and systems biology approaches. Methods: RNA-Seq data from peripheral blood mononuclear cells (PBMCs) of healthy persons, mild and severe 17 COVID-19 patients were analyzed to generate a gene expression matrix. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules in healthy samples as a reference set. For differential co-expression network analysis, module preservation and module-trait relationships approaches were used to identify key modules. Then, protein-protein interaction (PPI) networks, based on co-expressed hub genes, were constructed to identify hub genes/TFs with the highest information transfer (hub-high traffic genes) within candidate modules. Results: Based on differential co-expression network analysis, connectivity patterns and network density, 72% (15 of 21) of modules identified in healthy samples were altered by SARS-CoV-2 infection. Therefore, SARS-CoV-2 caused systemic perturbations in host biological gene networks. In functional enrichment analysis, among 15 non-preserved modules and two significant highly-correlated modules (identified by MTRs), 9 modules were directly related to the host immune response and COVID-19 immunopathogenesis. Intriguingly, systemic investigation of SARS-CoV-2 infection identified signaling pathways and key genes/proteins associated with COVID-19's main hallmarks, e.g., cytokine storm, respiratory distress syndrome (ARDS), acute lung injury (ALI), lymphopenia, coagulation disorders, thrombosis, and pregnancy complications, as well as comorbidities associated with COVID-19, e.g., asthma, diabetic complications, cardiovascular diseases (CVDs), liver disorders and acute kidney injury (AKI). Topological analysis with betweenness centrality (BC) identified 290 hub-high traffic genes, central in both co-expression and PPI networks. We also identified several transcriptional regulatory factors, including NFKB1, HIF1A, AHR, and TP53, with important immunoregulatory roles in SARS-CoV-2 infection. Moreover, several hub-high traffic genes, including IL6, IL1B, IL10, TNF, SOCS1, SOCS3, ICAM1, PTEN, RHOA, GDI2, SUMO1, CASP1, IRAK3, HSPA5, ADRB2, PRF1, GZMB, OASL, CCL5, HSP90AA1, HSPD1, IFNG, MAPK1, RAB5A, and TNFRSF1A had the highest rates of information transfer in 9 candidate modules and central roles in COVID-19 immunopathogenesis. Conclusion: This study provides comprehensive information on molecular mechanisms of SARS-CoV-2-host interactions and identifies several hub-high traffic genes as promising therapeutic targets for the COVID-19 pandemic.

Understanding Gene Expression and Transcriptome Profiling of COVID-19: An Initiative Towards the Mapping of Protective Immunity Genes Against SARS-CoV-2 Infection.

The COVID-19 pandemic has created an urgent situation throughout the globe. Therefore, it is necessary to identify the differentially expressed genes (DEGs) in COVID-19 patients to understand disease pathogenesis and the genetic factor(s) responsible for inter-individual variability. The DEGs will help understand the disease's potential underlying molecular mechanisms and genetic characteristics, including the regulatory genes associated with immune response elements and protective immunity. This study aimed to determine the DEGs in mild and severe COVID-19 patients versus healthy controls. The Agilent-085982 Arraystar human lncRNA V5 microarray GEO dataset (GSE164805 dataset) was used for this study. We used statistical tools to identify the DEGs. Our 15 human samples dataset was divided into three groups: mild, severe COVID-19 patients and healthy control volunteers. We compared our result with three other published gene expression studies of COVID-19 patients. Along with significant DEGs, we developed an interactome map, a protein-protein interaction (PPI) pattern, a cluster analysis of the PPI network, and pathway enrichment analysis. We also performed the same analyses with the top-ranked genes from the three other COVID-19 gene expression studies. We also identified differentially expressed lncRNA genes and constructed protein-coding DEG-lncRNA co-expression networks. We attempted to identify the regulatory genes related to immune response elements and protective immunity. We prioritized the most significant 29 protein-coding DEGs. Our analyses showed that several DEGs were involved in forming interactome maps, PPI networks, and cluster formation, similar to the results obtained using data from the protein-coding genes from other investigations. Interestingly we found six lncRNAs (TALAM1, DLEU2, and UICLM CASC18, SNHG20, and GNAS) involved in the protein-coding DEG-lncRNA network; which might be served as potential biomarkers for COVID-19 patients. We also identified three regulatory genes from our study and 44 regulatory genes from the other investigations related to immune response elements and protective immunity. We were able to map the regulatory genes associated with immune elements and identify the virogenomic responses involved in protective immunity against SARS-CoV-2 infection during COVID-19 development.

An integrative bioinformatics analysis for identifying hub genes associated with infection of lung samples in patients infected with SARS-CoV-2.

BACKGROUND: At the end of 2019, the world witnessed the emergence and ravages of a viral infection induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Also known as the coronavirus disease 2019 (COVID-19), it has been identified as a public health emergency of international concern (PHEIC) by the World Health Organization (WHO) because of its severity. METHODS: The gene data of 51 samples were extracted from the GSE150316 and GSE147507 data set and then processed by means of the programming language R, through which the differentially expressed genes (DEGs) that meet the standards were screened. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on the selected DEGs to understand the functions and approaches of DEGs. The online tool STRING was employed to construct a protein-protein interaction (PPI) network of DEGs and, in turn, to identify hub genes. RESULTS: A total of 52 intersection genes were obtained through DEG identification. Through the GO analysis, we realized that the biological processes (BPs) that have the deepest impact on the human body after SARS-CoV-2 infection are various immune responses. By using STRING to construct a PPI network, 10 hub genes were identified, including IFIH1, DDX58, ISG15, EGR1, OASL, SAMD9, SAMD9L, XAF1, IFITM1, and TNFSF10. CONCLUSION: The results of this study will hopefully provide guidance for future studies on the pathophysiological mechanism of SARS-CoV-2 infection.

Transcriptomic Profiles Reveal Downregulation of Low-Density Lipoprotein Particle Receptor Pathway Activity in Patients Surviving Severe COVID-19.

To assess the biology of the lethal endpoint in patients with SARS-CoV-2 infection, we compared the transcriptional response to the virus in patients who survived or died during severe COVID-19. We applied gene expression profiling to generate transcriptional signatures for peripheral blood mononuclear cells (PBMCs) from patients with SARS-CoV-2 infection at the time when they were placed in the Intensive Care Unit of the Pavlov First State Medical University of St. Petersburg (Russia). Three different bioinformatics approaches to RNA-seq analysis identified a downregulation of three common pathways in survivors compared with nonsurvivors among patients with severe COVID-19, namely, low-density lipoprotein (LDL) particle receptor activity (GO:0005041), important for maintaining cholesterol homeostasis, leukocyte differentiation (GO:0002521), and cargo receptor activity (GO:0038024). Specifically, PBMCs from surviving patients were characterized by reduced expression of PPARG, CD36, STAB1, ITGAV, and ANXA2. Taken together, our findings suggest that LDL particle receptor pathway activity in patients with COVID-19 infection is associated with poor disease prognosis.

HumanNet v3: an improved database of human gene networks for disease research.

Network medicine has proven useful for dissecting genetic organization of complex human diseases. We have previously published HumanNet, an integrated network of human genes for disease studies. Since the release of the last version of HumanNet, many large-scale protein-protein interaction datasets have accumulated in public depositories. Additionally, the numbers of research papers and functional annotations for gene-phenotype associations have increased significantly. Therefore, updating HumanNet is a timely task for further improvement of network-based research into diseases. Here, we present HumanNet v3 (https://www.inetbio.org/humannet/, covering 99.8% of human protein coding genes) constructed by means of the expanded data with improved network inference algorithms. HumanNet v3 supports a three-tier model: HumanNet-PI (a protein-protein physical interaction network), HumanNet-FN (a functional gene network), and HumanNet-XC (a functional network extended by co-citation). Users can select a suitable tier of HumanNet for their study purpose. We showed that on disease gene predictions, HumanNet v3 outperforms both the previous HumanNet version and other integrated human gene networks. Furthermore, we demonstrated that HumanNet provides a feasible approach for selecting host genes likely to be associated with COVID-19.

New insights into genetic characteristics between multiple myeloma and COVID-19: An integrative bioinformatics analysis of gene expression omnibus microarray and the cancer genome atlas data.

BACKGROUND: Multiple myeloma (MM) is a hematological malignancy. Coronavirus disease 2019 (COVID-19) infection correlates with MM features. This study aimed to identify MM prognostic biomarkers with potential association with COVID-19. METHODS: Differentially expressed genes (DEGs) in five MM data sets (GSE47552, GSE16558, GSE13591, GSE6477, and GSE39754) with the same expression trends were screened out. Functional enrichment analysis and the protein-protein interaction network were performed for all DEGs. Prognosis-associated DEGs were screened using the stepwise Cox regression analysis in the cancer genome atlas (TCGA) MMRF-CoMMpass cohort and the GSE24080 data set. Prognosis-associated DEGs associated with COVID-19 infection in the GSE164805 data set were also identified. RESULTS: A total of 98 DEGs with the same expression trends in five data sets were identified, and 83 DEGs were included in the protein-protein interaction network. Cox regression analysis identified 16 DEGs were associated with MM prognosis in the TCGA cohort, and only the cytochrome c oxidase subunit 6C (COX6C) gene (HR = 1.717, 95% CI 1.231-2.428, p = .002) and the nucleotide-binding oligomerization domain containing 2 (NOD2) gene (HR = 0.882, 95% CI 0.798-0.975, p = .014) were independent factors related to MM prognosis in the GSE24080 data set. Both of them were downregulated in patients with mild COVID-19 infection compared with controls but were upregulated in patients with severe COVID-19 compared with patients with mild illness. CONCLUSIONS: The NOD2 and COX6C genes might be used as prognostic biomarkers in MM. The two genes might be associated with the development of COVID-19 infection.

Network Pharmacology Reveals That Resveratrol Can Alleviate COVID-19-Related Hyperinflammation.

Hyperinflammation is related to the development of COVID-19. Resveratrol is considered an anti-inflammatory and antiviral agent. Herein, we used a network pharmacological approach and bioinformatic gene analysis to explore the pharmacological mechanism of Resveratrol in COVID-19 therapy. Potential targets of Resveratrol were obtained from public databases. SARS-CoV-2 differentially expressed genes (DEGs) were screened out via bioinformatic analysis Gene Expression Omnibus (GEO) datasets GSE147507, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis; then, protein-protein interaction network was constructed. The common targets, GO terms, and KEGG pathways of Resveratrol targets and SARS-CoV-2 DEGs were confirmed. KEGG Mapper queried the location of common targets in the key pathways. A notable overlap of the GO terms and KEGG pathways between Resveratrol targets and SARS-CoV-2 DEGs was revealed. The shared targets between Resveratrol targets and SARS-CoV-2 mainly involved the IL-17 signaling pathway, NF-kappa B signaling pathway, and TNF signaling pathway. Our study uncovered that Resveratrol is a promising therapeutic candidate for COVID-19 and we also revealed the probable key targets and pathways involved. Ultimately, we bring forward new insights and encourage more studies on Resveratol to benefit COVID-19 patients.